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Abstract 

We study instanton effects on inclusive semileptonic b — ► u decay using the 
heavy quark effective field theory and the operator product expansion. The 
effect contributes not only to the hadronic matrix element but also to the 
coefficient functions in the operator product expansion. We find that the 
coefficient function is singular near the boundaries of the phase space. In 
order to use perturbative QCD reliably, it is necessary to introduce smearing 
near the boundary. However the instanton contribution to b — > c decay seems 
negligibly small. 
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I. INTRODUCTION 



The Cabibbo-Kobayashi-Maskawa matrix element V u b may be extracted from the shape 
of electron energy spectrum in inclusive semileptonic B decay. In order to avoid large b — > c 
decay background, we have to look at the endpoint region of (9(300MeV). However this 
is precisely the regime in which theoretical prediction is least understood due to various 
low-energy QCD effects. Therefore experimental data have been compared mainly with 
phenomenological models such as ACCMM or ISGW O. 

Recently there has been important theoretical progress in inclusive semileptonic B decay. 
The first model-independent QCD approach has been formulated using the heavy quark 
effective field theoy(HQEFT) and the operator product expansion(OPE). Chay, Georgi and 
Grinstein |3j have shown that the weak decay of B mesons can be described by a systematic 
expansion in inverse powers of m^, in which the leading contribution reproduces the parton- 
model result (See also Ref. @.). Using the OPE, Chay et.al. have shown that next-order 
corrections are of order (A/m?,) 2 , where A is the QCD scale, hence expected to be small. 
Mannel 0, Manohar and Wise || have analyzed these corrections and concluded that 
these corrections are singular at the boundaries of the Dalitz plot, which necessitates some 
smearing prescription before comparing with experimental data. 

There are also radiative QCD corrections to the decay rate J7|. In the leading log approx- 
imation, double logarithms are summed to provide the Sudakov factor. In b — > u decay, this 
changes significantly the shape of the electron spectrum near the endpoint. Falk, Jenkins, 
Manohar and Wise || have considered the validity of this resummation near the endpoint 
of electron energy. By studying the pattern of the summation of subleading logarithms, 
they have reached a conclusion that a smearing is necessary for the validity of the leading 
log approximation. Bigi et.al. || have reached essentially the same conclusion but with a 
careful study of running a s {k 2 ). 

The idea of using the OPE in the approach of Ref. || is that it factorizes the short- 
distance and the long-distance physics in terms of the product of the Wilson coefficient 



2 



functions and the matrix elements of local operators. So far the main focus of recent dis- 
cussions has been on the systematic analysis of nonperturbative hadronic matrix elements. 
However the coefficient functions may be affected by nonperturbative effects as well. One 
of the examples is an instanton contribution. 

In this paper we investigate the effect of instantons on the coefficient functions in B decay. 
Due to asymptotic freedom a s (Q 2 ) ~ 1/ InQ 2 as Q 2 — > oo, we expect that the strength of a 
single instanton is of order exp(— 2-n /a s (Q 2 )). At large momenta this is negligible compared 
to perturbative QCD corrections. However the characteristic scale for a s is not but 
the momentum transfer to final hadrons. Therefore near the endpoint of electron energy 
spectrum where a s becomes large, it is likely that instanton effects on the decay rate may 
grow sizably. For large momentum transfer Q 2 to final hadrons, we expect that small 
instantons of size p ^ 1/\Q\ are most relevant. We study their effect by calculating one 
instanton (+ anti-instanton) corrections to the inclusive semileptonic b decay rate. We find 
that the effect is negligible in most of the the phase space except near the boundaries of the 
Dalitz plot. At the boundaries we find that the instanton contribution becomes singular. 
Therefore the decay rate has to be smeared in such a way that the instanton contribution 
remains small compared to the parton-model result. 

In this estimate we use the dilute gas approximation since the main contribution comes 
from small instantons. Instantons contribute not only to the hadronic matrix elements but 
also to the coefficient functions. We emphasize that we have neither attempted to calculate 
QCD radiative corrections in the realistic instanton background nor other nonperturbative 
effects such as multi-instantons or renormalons. It is possible that any of these give bigger 
contributions than a single instanton contribution. Even in this case our estimate may still 
characterize a typical size of nonperturbative effects in the coefficient functions. After all 
the dilute instantons are the only known calculable nonperturbative effect in QCD. 

This paper is organized as follows. In Section 2 we discuss the kinematics relevant to the 
inclusive semileptonic B decay. In Section 3 we recapitulate the OPE and summarize the 
general structure of local operators and their coefficient functions. In Section 4 we compute 
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the instanton contribution to the coefficient functions in the forward Compton scattering 
amplitude. In Section 5 we discuss numerical estimate of the results of Section 4. We 
conclude in Section 6 and also comment on the nonperturbative effects in other cases. 



II. KINEMATICS 

The semileptonic b — > u decay is described by the weak Hamiltonian 

G F 

H w = ^=V^ti7 M (l + 75)6 ■ ej^l + 75K 

= ^=V ub JMJ+,, (1) 

where J£ = ^(1 + 75)6 and Jg = P e 7^(l + 7 5 )e. For b — > c decay the Hamiltonian is obtained 
from Eq. ([!]) by replacing u and V ub by c and V cb respectively. The inclusive decay rate is 
related to the hadronic tensor defined by 

Wr = (2K)'Y,P(PB-q-Px)(B\J!*\X)(X\JZ\B), (2) 
x 

where X denotes all possible hadronic final states. Let k£ be the electron momentum, k£ 
be the antineutrino momentum, and p B = Mbv^ be the momentum of the B meson. Then 
the hadronic tensor depends on the timelike vector g M = k£ + hp. In this paper we do 
not distinguish the B meson mass Mb and the b quark mass m b . This difference may be 
incorporated into the kinematics as discussed in Refs. Hf§. 
The differential decay rate of the B meson is given by 

dT = -^\V ub \ 2 W^L llu d(PS), (3) 
4m b 

where L^ v is the leptonic tensor given by 

IT = 8(A£A£ + k v e k>" 9 - gTK ■ h - ie» v aP k«kl). (4) 

We can express the phase space integral in terms of the three independent kinematic variables 
y = 2E e /m b , q 2 = q 2 jm\ and v ■ q = v ■ q/m b as 



In the complex v ■ q plane with q 2 fixed, W^ v is related to the discontinuity of the forward 
Compton scattering amplitude T^ u across a physical cut as 

^ = 2Imr, (6) 

where 

T^(q,v) = -ij d'xe^ (B\T{ J^r u (0)}\B). (7) 

To evaluate the decay rate it is necessary to examine the analytic structure of T^(g, v ) 
|3|||. In the complex v ■ q plane, a "physical cut" relevant to the decay is located on the 
real axis between \fq^ < v ■ q < (1 + q 2 — p)/2 where p = m 2 /m 2 is the mass ratio- 
squared of the u quark to the b quark. Discontinuity of T^ u across this physical cut yields 
W^ u for semileptonic B decay. There are also other cuts for v ■ q < —\fq I and v ■ q > 
|((2 + ^/p) 2 — q 2 — 1) but they correspond to other physical processes. In Fig.l we show the 
analytic structure of T^ u in v ■ q plane. As v ■ q approaches (1 + q 2 — p 2 )/2 the final hadron 
reaches a resonance region. There we expect significant nonperturbative QCD effects. Chay 
et.al. have observed that one can calculate the decay rate even in this limit perturbatively 
by choosing the integration contour away from the resonance region. This should provide a 
good approach to b — > c inclusive decay. Since m 2 ^> A 2 , the other cut on the right-hand 
side in Fig.l always stays away from the physical cut. 

However this is no longer true in b — ■> u decay. In this case we cannot choose a contour 
away from the resonance region without enclosing some of the unphysical cut. When q 2 — ► 1 
and p — > 0, the physical cut and the right-hand side cut in Fig.l pinch together. This is 
in fact the origin of large logarithms of lepton energy found in Ref. 0. Since the charm 
threshold lies at y — 1 — p ~ 0.9, we have to deal with this problem. For the moment we 
consider q 2 m 2 but (m^v — q) 2 ^> A 2 in order to calculate the decay distribution reliably 
using perturbative QCD. In the end we are interested in how far the result may be extended 
to the endpoint region. 
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The total decay rate is given by 

r = l y *l y d » /o rfg 7| + £^ ■ g (8) 

Note that the interval of v ■ q integration in Eq. (||) does not coincide with the physical cut. 
From the kinematics the minimum of v ■ q is at point P in Fig.l where v ■ q = y/2 + q 2 /2y. 
In Eq. (H) W^L^y is evaluated from the discontinuity of T^L^ along the contour C in 
Fig.l. Using Cauchy's theorem this is related to the contour integral along the contour C in 
the complex v ■ q plane. The integral along C is reliably evaluated using perturbative QCD 
as long as (m&i; — q) 2 ^> A 2 . Then the double differential decay rate can be expressed as 

d 2 T G 2 F m\ |T „ „ f , A T^L 



rfg 2 rfy 2 8 tt 4 



u6 1 



/ dv ■ q (9) 

Jc mi 



The contour C is shown in Fig.l. Care must be taken in choosing the contour C in Eq. @. 
For example if T^ v contains poles only, any contour enclosing the poles gives the right 
answer 0,0. On the other hand if T^ u has cuts, the contour has to be chosen in such a 
way to cover the correct v ■ q integration interval determined by kinematics. We therefore 
have chosen the contour C in Eq. (|9|) as a circle of radius z = (y — q 2 )(l — y)/y centered at 
v-q = (l+g 2 )/2. The contribution of the contour integral near the point P, v-q = y /2+q 2 /2y 
is negligibly small as long as 2 > A/m^. 

Near the boundaries of the Dalitz plot q 2 — > y or y — > 1, z shrinks to zero and we are in 
the resonance region. This is where nonperturbative effects dominate. This limits the extent 
of our theoretical prediction of the decay rate to which the experimental data can be com- 
pared. Therefore it is necessary to introduce a smearing in order to use perturbation theory 
reliably JK|. The size of the smearing is determined by requiring that the nonperturbative 
effects be smaller than the parton- model result. 



III. OPERATOR PRODUCT EXPANSION 



For processes involving 6-quark decay, it is appropriate to use the HQEFT. In the 
HQEFT, the full QCD b field is expressed in terms of b v for b quark velocity v. The b v 
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field is defined by 



b v = l±L e im b vx b + . . . ^ ( 1Q j 



where the ellipses denote terms suppressed by powers of l/m&. The equation of motion for 
b v is fb v = b v . We can apply the techniques of the OPE to expand T^ v in terms of matrix 
elements of local operators involving b v fields 



= -il dW*-*(£|T{J£(a;)tj£(0)}|B) 

= Y.C» v {v,q)(B\0^\B). (11) 

n,v 

Here are local operators involving b v b v bilinears. The coefficient functions C£ u depend 
explicitly on q and v because of the v dependence of O^ 1 '. All large momenta are contained 
in the coefficient functions while the matrix elements describe low-energy physics. 

In this scheme T^ u is expanded as a double series in powers of a s and l/m&. To leading 
order in a s and l/m b 



T^ = -i j d 4 xe^ mbV -^ x (B\T{b v (x)^(l + l5 )u(x)u(0)Y(l + -y s )b v (0)}\B) 

= -i [d 4 xe^-^ x (B\T{b v (x)Y(l+l 5 )S (x)Y(l + l 5 )b v (0)}\B), (12) 



where Sq(x) denotes a free u quark propagator. In momentum space Eq. ( |T2| ) is proportional 
to 

^(§TW^ 1 + ^ < 13) 

where Q = m^v — q is the momentum transferred to the u quark and is a small residual 
momentum of the b v field of order A. 

When A 2 < Q 2 C ml, it is sufficient to keep the terms at leading order in l/m& only 
and expand Eq. ([13]) in powers of k/\Q\. This generates coefficient functions and local 
operators in which k is replaced by the derivatives acting on the b v fields. The leading term 
independent of k gives the parton-model result. Contracting the leading term with L^ u we 
find 



ITV = 64^^(516^(1 + l5 )b v \B) 

= l2Sm b K - Q £-\ (14) 

where T^ u is T^ u at leading order in k/\Q\. The second line in Eq. (|14]) follows from the 
normalization (B\b v b v \B) = 2ra&. 

The next-order correction to Eq. ([14]) is suppressed by at least (A/Q) 2 0. Corrections 
to from higher dimensional operators involving b v b v bilinears have been studied system- 
atically Furthermore Neubert [|11]] has resummed these corrections to get a "shape 
function". He has observed that the shape function is universal independent of the final 
quark flavor. 

Usually the coefficient functions are calculated perturbatively. In addition there could 
be nonperturbative correction to the coefficient functions as well ||12|| . In the next section 
we calculate how instantons affect the coefficient function of the leading term, Eq. (|14|). 
in the OPE. The instantons contribute not only to the matrix elements but also to the 
coefficient functions. If there is an infrared divergence due to large instantons we attribute 
it to the contribution to the matrix elements. For infrared-finite parts we interpret them as 
the correction to the coefficient functions. 



IV. INSTANTON CONTRIBUTION 

We now compute the contribution of instantons to the coefficient functions. In estimating 
the contribution we start from the Euclidean region where Q 2 is large enough to use the 
OPE reliably. We expect that the main contribution comes from small instantons of size 
p ^ 1/121' hence we use the dilute gas approximation in what follows. More specifically we 
calculate the instanton correction to the decay rate at leading order in both a s and k in the 
OPE. This is the instanton correction to the parton model result Eq. (|14T). 

In the background of an instanton (+ anti-instanton) of size p and instanton orientation 
U located at the origin, the Euclidean fermion propagator may be expanded in small fermion 
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masses as ||13 



S±{x, y; p±; U±) = -—if>o(x)if>l(y) + S^\x, y\ p±; U ± 
+ mj d A wS { 2 ) (x,w;p±]U±)S { ± ) (w,y;p±;U±) 



(15) 



where ± denotes instanton, anti-instanton. ?/> is the fermion zero mode eigenfunction and 
S±^ = J2e>o E±{x)^ E ±{y) is the Green's function of fermion nonzero modes. 

In evaluating the forward Compton scattering amplitude T^ v , Eq. (|7|), we should use the 
propagator in Eq. (|T5|) instead of Sq. T^ v can be written as 



1 J inst.5 



(16) 



where the first term is the parton-model amplitude. The second term is the amplitude due 
to instantons of all orientation U, position z and size p. A schematic configuration of an 
instanton or anti-instanton is shown in Fig. 2. After averaging over instanton orientations 

T inst is S iven b Y 



inst. 



rf 4 Ae , S .A [d*z a dp a D(p a )(B\b v (x)j»(l + l5 ){S a (X,Y;p a ) - S (A)}Y(l + -ft)b v {y) \B), 

a=± J 



(17) 



where D(p) is the instanton density, A = x — y, X = x — z and Y = y — z. In Eq. (|T] 
S±(X, Y; p±) is the fermion propagator averaged over instanton (anti-instanton) orientations 
centered at z. Using the MS scheme with rif flavors of light fermions, D(p) is given by |14[ 



D{p) = KA b (pA) b+ 3 In 



p 2 A 2 



(18) 



where 



X 



7T 



exp[- - a(l) + 2(n/ - l)a(- 



(19) 
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in which the (3 function at two loops and the running mass at one loop are used and rhi 
are the renormalization-invariant quark masses. In Eq. ( |i~9|) ot(l) = 0.443307 and a(l/2) = 
0.145873. From now on we replace the logarithmic term in D(p) by its value for p = 1/|Q|. 
Corrections to this replacement are negligible since they are logarithmically suppressed. 

Inserting Eq. (|T^) to Eq. flT7|), the leading contribution comes from the mass-independent 
part S± due to the chiral structure of the left-handed weak currents. The chirality- 
conserving part of S±^ in singular gauge |13| is written as 



E sS-\x, Y; p) = Jr(x 2 r 2 + P 2 x ■ y)[x 2 y 2 (x 2 + P 2 )(y 2 + P 2 )]^ 



a=± 

4vr 2 A 2 



(X 2 [X 2 F 2 (X 2 + p 2 f{Y 2 + p 2 )}-i + (X <-> Y)) 



p2 (X [X 2 F 2 (X 2 + p 2 )\Y 2 + p 2 )]-5 _ (x <-> Y)). (20) 



4tt 2 

From Eq. (p0|) S±(X,Y; p) — S (A) can be written as 

E P) - 5 (A)) « -|^^(X 2 - F 2 ) 2 [X 2 F 2 (X 2 + p 2 )(F 2 + p 2 )]- 1 



+ ^l(x^7 + V^T? ) [X2Y2{X2 + p2){Y ' 2 + p2)X 
4tt 2 vx 2 + p 2 r 2 + p 2 



' [X 2 F 2 (X 2 + p 2 )(F 2 + p 2 )]-5, (21) 



where Sq(A) = —A , /2tt 2 A 4: is the free quark propagator. In deriving Eq. ([H]) we move 
— p 2 (X 2 + Y 2 )/2 proportional-part of the first line to the second line in Eq. ( p0|) using 
p 2 X • K = p 2 ((X + Y) 2 — X 2 — y 2 )/2. We also approximate the remaining terms to get 
the first line in Eq. (|21|) so that exact analytic calculation is possible while leaving small 
instanton contribution essentially unaltered. Plugging Eq. ( pip into Eq. (fL7l), we integrate 
over the instanton center z, the instanton size p, and finally make a Fourier transform over 
A. Note that the integral over z is convergent as \z\ — > oo. 

We notice that the integration over p is convergent for small instantons. On the other 
hand large instanton part is divergent. However since the integrand is analytic for large p 
(See Appendix.), there are only a finite number of divergent terms when T^ st is expanded 
in 1/p. We interpret these infrared divergent terms as the instanton contribution to the 
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matrix elements of operators in the OPE |15|. The remaining terms are infrared convergent 



and are interpreted that they contribute to the coefficient functions. In order to calculate 
finite terms any convenient regularization prescription will do. To this end we analytically 
continue the exponent of p in the instanton density so that D(p) oc p M ~ 4 |16|] and the 
spacetime dimensions to d = 4 + 2e. In the end we let M — > 11 for rif = 3 and e — > to get 
the final answer. 

After some algebra which we show in the Appendix, the contribution to T^L^ from 
the first line of Eq. (^) consists of two pieces. The contribution of the first piece is written 

as 



Kjdpp 4 D(p) J 



d 4 w 



w 2 (w 2 + p 2 ) 



x 2 6 ^Lk^B\b vl ,{l + l5 )b v \B). (22) 

The integration over p is divergent as we take the physical value M = 11. This is interpreted 
as the large-instanton contribution to the matrix element b v ^^b v . 

The contribution from the second piece is regular as M — > 11, e — ► and is given by 

2 10 2 2/~n M 2 \ 6 A Q 2 \ 10/9 



x 2 6 ^—±k^(B\b vlll (l + lh )b v \B). (23) 



,6 fee • Q 

Q 2 

The contribution from the second and the third lines in Eq.(|2l|) to T^L^ vanishes, as 
explained in the Appendix. Therefore Eq. (|23| ) is the overall finite contribution. 

It is now straightforward to evaluate the instanton contribution to the differential decay 
rate. So far we have worked in Euclidean spacetime. We now make a naive analytic contin- 
uation to Minkowski spacetime with timelike Q 2 = (m^v — q) 2 to evaluate the differential 
decay rate. Then T^L^ is written as 

-wj _ 5^ 2 (vq-l)(y-q 2 ) 



ln^(^-g-^-)_ 
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10/9 



(24) 



where 



2 n 7r 2 / A \ 12 

A -iw r ^ K ■ (25) 



Eq. ([24] ) has a branch cut emanating from the resonance point v ■ q = 1 + q 2 /2, which we 
have shown in Fig.l. 

To get the double differential decay rate d 2 T/dq 2 dy, we integrate Eq. (]24]) over v-q along 
C in Fig.l as described in Section 2. The differential decay rate is written as 

1 rfr inst 15 A fj a (y-q 2 )(vq-l) 

S2> 



-A \ dv ■ q 

■ Jc 



T dq 2 dy 32vr 7c (« • q - i±2l)7 



x 



\n—^-(v-q 



2 

,2 i , -2 -,10/9 



-2mg A 1 + g 2 



(26) 



where T = G 2 ?] V u i,\ 2 ml /192tt 3 . The exponent of the logarithm is almost unity and we 
replace 10/9 by 1. We expect that it does not change the result much. We finally get 

1 ci 2 r inst , 5 g 2 -(l-y)(y-g 2 ) 



T rfg 2 % ' (1 - y) 6 (y - q 2 f 
The instanton effect is suppressed at y ~ 0. However it become singular near the boundaries 
of the phase space at y ~ 1 and y ^ q 2 ■ At these boundaries the final quark approaches the 
resonance region. As there are large instanton contributions near the resonance point, it is 
necessary to introduce a smearing to define sensible decay rate in perturbative QCD. We 
discuss this in detail in the next section. 



V. NUMERICAL ANALYSIS 

We have calculated the instanton effect on the inclusive semileptonic b decay. As we 
have observed in the previous section the effect is singular at the boundaries of the phase 
space. In order to use perturbative expansions reliably a smearing prescription has to be 
introduced. 

As a simple prescription we first consider the smeared single differential decay rate 

<^>^/V % -f-^ (28) 
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in which we restrict the phase space so that the singular region y = q 2 is avoided by 
5. The size of smearing 5 is determined by the requirement that the smeared instanton 
contribution to the single differential decay rate be smaller than that in the parton model 
dTo/dy = 2y 2 {3 - 2y)T . Eq.(f§) is given by 

/ Gainst dT A 1 



x dy v ^ y 24(l- 2 /)6(3-2y) 

x [9-4y-24(|) 3 + 15(|) 4 + 4(^)]. (29) 

The ratio R(y, 5) = (dT inst /dy)s/(dT /dy) is plotted in Fig.3 for three different values of 
5 = 0.15, 0.17, 0.19 with A = 400 MeV and mj = 5 GeV. The numerical value of A is given 
by 

V m b J \m b J 



in which we set nt = 3 and the renormalization- invariant quark masses as |17 



m u = 8.2 ± 1.5MeV, 
m d = 14.4 ± 1.5MeV, 

m s = 288 ± 48MeV. (31) 

We see that for these choices of 5, R(y, S) grows rapidly for y > 0.84. From Fig.3 we see 
that a smearing of 5 ~ 0.16 ~ 0.20 at the boundary y = q 2 of the phas space is needed to 
extract a reliable electron spectrum from the single differential decay rate. 

We have also examined the behavior of R(y, 5) as we increase 5. The position at which 
R(y, 5) grows like a brick wall does not shift much from y ~ 0.84. Because of this singular 
behavior we need another smearing near y — 1. For simplicity we cut off the region near 
y = 1 by the same smearing width 5 as in Eq. (|2T|). 

Instanton contribution to the total decay rate after such a smearing prescription is given 

by 

OW, = Cdy [ V dq 2 9(y - q 2 - 5)9(1 -y- 8) dH iuM 



o Jo dq 2 dy 
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r ^7 4( 6 - 535 + 1985 2 - 4205 3 + 6125 4 



A 1 

24 ¥' 

3545 5 + 165 6 - 45 7 - 5 9 + 1805 5 ln<5). (32) 



The ratio R(5) = (T) 5 /T is plotted in Fig.4 for three different values of A = 350,400,450 
MeV and for nit, = 5 GeV. We see that the ratio R(S) also rises sharply like a brick wall as 5 
decreases. When 5 is large, say 5 > 0.15, R(5) is insensitive to the choice of A/m^. However, 
for small 5, say 5 ~ 0.12, R(S) is sensitive to the value of A/rrib. If we require R(S) be less 
than 20%, the size of the smearing is roughly at least 0.12 ~ 0.15 for our choices of A. 
Our results Eqs. (P^), (|3"2"D are sensitive to the precise value of A. This is what we expect 



from the instanton effect. Nevertheless, as emphasized in Introduction, our results may still 
represent a typical size of nonperturbative effects. 



VI. CONCLUSION 

In this paper we have studied instanton effects on the decay rate of the inclusive semilep- 
tonic b — > u decay. In particular we have estimated the instanton correction to the coefficient 
function at leading order. We have found that the correction is singular at the boundaries 
of the phase space. Because of the large instanton contribution to the coefficient function, 
we have found that it is necessary to introduce a smearing near the boundaries of the phase 
space. For the single differential decay rate a smearing of size 5 ~ 0.16 ~ 0.2 is needed 
taking into account of the correlation of the smearing sizes at the two boundaries y = 1 and 
y = q 2 . For the total decay rate the smearing size is 5 ~ 0.12 ~ 0.15. Note that dT^/dy 
is appreciable for all y while dT inst /dy is negligible except at the boundaries of the Dalitz 
plot. Therefore the smearing size for the total decay rate is smaller than that for the single 
differential decay rate. 

Precise theoretical prediction in the endpoint region is important for the extraction of V u b. 
Singular nature of various corrections indicates we need some smearing before comparing 
with data. Manohar and Wise have studied the smearing necessary for the effects of 
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matrix elements of higher dimensional operators. They have concluded that the size of 



smearing should be e ~ 0.2. Neubert |TT| has suggested an idea to measure the universal 
shape function from b — > c decay and apply it to b — > u decay. He has proposed this as a 
promising way to a precise extraction of V u b- 

Falk et.al. || have studied the smearing for the validity of the leading log approximation 
and have found the smearing size should be roughly 0.1 or less. The instanton effect is 
another source of corrections to B decays. While it is negligible for b — > c decay the 
effect is significant in b — > u decay. Our analysis shows that the smearing has to be about 
5 ~ 0.12 ~ 0.16. This is comparable to the difference in the endpoints of the b —>■ c and 
b — > u decays. Therefore it is perhaps difficult to eliminate a model dependence in extracting 
V u b from inclusive semileptonic b — > u decay. 

We can similarly consider the instanton effect in inclusive semileptonic c decay. As we 
can infer from the quark mass dependence in Eq. (|25"D we expect the effect is much larger 
in this case because m c <C mj. Therefore the decay rate near the endpoint is strongly 
model-dependent due to large uncertainties from radiative corrections and nonperturbative 
effects H including instantons. We may also consider the instanton effect on inclusive b — > c 
decay. Because Q 2 > m 2 . the instanton effect is finite even at the boundaries of the phase 
space. The numerical value is negligibly small since it is proportional to high powers (A/m c ) n 
where n may be determined from the instanton density Eq. (|18D for rif = 4. 
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APPENDIX: INSTANTON INTEGRATIONS 



There are three integrations [See Eq. (21).] we have to perform. The integrals of Eq. (|2T 
over p, z and Fourier transform of A after a suitable change of variables are expressed as 

/ O 2 \ 10 / 9 r°° r n 2 

Klln^-) [ dpp M ~ A d d td d u P 



x 



A 2 J Jo rr J 4vr 2 

(33) 



ide iQt \f p 2 p 2 



d$ t 4 J\u 2 (u 2 + p 2 ) u 2 ((t - u) 2 + p 2 ) 
id e iQ t \ 2p 2 

d $ t2 ' ^u 2 {u 2 + p 2 )(t - u) 2 ((t - u) 2 + p 2 ) 

id e iQ - f 



(34) 

(35) 

' d $ ^t\t 2 + P 2 yJ ^u 2 (u 2 + P 2 ) 

where 8/80 = ^y^d/dQ' 1 . At small p, the p integral is convergent while at large p it is 
divergent. However for large p the integrand may be expanded in powers of 1/p. We can 
see clearly that there are only a finite number of divergent terms at the physical value 
M = 11 for nf = 3. These divergent terms are interpreted as the instanton contribution to 
the matrix elements of operators in the OPE [|T5|]. Since the remaining terms are finite we 
may evaluate them by any convenient methods. We have analytically continued M and the 
spacetime dimensions d = 4 + 2e. We set M = 11 and d = 4 in the end. 

The first term in Eq. (|33| ) is divergent and we interpret it as a contribution to the hadronic 
matrix elements |15| . In evaluating the second term in Eq. (|33"D we first integrate over u and 
then over p. Finally we make a Fourier transform with respect to t. The result turns out to 
be regular as e — > 0. Setting e = the second term in Eq. (^) is given by 

w , 9 M + 1 9 /M+ i\r(- 1+A/x 



10 / O 2 \ 10 /9 

x^fl^) • (36) 



M + 3 V 2 Jr(-¥=Z) 

QM+3 V m A?) 

In getting Eq. (|36|) we have analytically continued M in such a way that the integral becomes 
convergent. However the result can be extended to M = 11 as it is convergent for any positive 
M > 3. The final result is 
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(M + 3)(M-l) r I— J I 1 " W ' (37) 



Putting M = 11 yields 



4.16xl0 6 K^(ln^j . (38) 



Eq. ( |3~4] ) is regular as e — > 0. Setting e = and after some calculation it is written as 

xF(H±I,J£±V + 2;l-|), (39) 

where F = 2-^1 is the hypergeometric function. Here x and y are Feynman parameters 
introduced in integrating over the instanton position. Eq. ( |3"5| ) is also regular as M — > 11, 
hence we can set M = 11. We can expand the hypergeometric function as a series of (y/x) n 
for n > 0. The y integral is of the form 

^{l-y)- 1 ' 2 , (40) 



Jo 



which is finite for all n > 0. On the other hand the x integral is of the form 

( l dxx- M / 2 ~ n -\l-xY l /\ (41) 

JO 

which is zero for odd M and infinity for even M. Therefore for M = 11 the integral vanishes. 

Similarly Eq. (|35|) is regular as e — > and M — > 11. After some calculation the result 
with e = 0,M = llis written as 



)12 r(7)r 2 (6)r(5) %<£ ^_ g 2 y°/9 



T(12) q 14 v a 2 / 
x T dxdyx- u ' 2 {l - x)- 1 ^- 172 !! - Z/) 1/2 





x F(5,6;12;l- -). (42) 

Again the hypergeometric function can be expanded as a series of (y/x) n for n > 0. The y 
integral is finite but the x integral vanishes for the same reason as in Eq. (f4"T|). Therefore it 



turns out that Eqs. (|34j) , (j35|) vanish and the only finite contribution is Eq. (|38f) . 
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FIGURES 

FIG. 1. Anaytic structure of T^ v in the complex v • q plane. The point P is the minimum of 
v • q at y/2 + q 2 /2y for b — > u decay. The contour integral along C' is related to the contour integral 
along C. 

FIG. 2. Schematic diagram of an instanton configuration with size p located at contributing 
to the forward Compton amplitude. 

FIG. 3. R(y,5) for 5= 0.15, 0.17, 0.19 with A = 400 MeV and m b =5 GeV. 

FIG. 4. R(6) for A= 350, 400, 450 MeV with m b = 5 GeV. 
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